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Abstract 

We propose a scanning tunneling microscopy Aharonov-Bohm (AB) interferometer on the surface 
of a topological insulator (TI) to probe the crossover from analogous weak antilocalization (WAL) 
to weak localization (WL) phenomenon via the AB oscillations in spin-resolved local density of 
states (LDOS). Based on our analytical and numerical results, we show that with increasing the 
energy gap of TI surface states, the <&o/2=hc/2e periodic AB oscillations in spin-resolved LDOS 
gradually transit into the <3?o periodic oscillations. 
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Topological insulators (TI) have attracted substantial interest in the modern condensed 

nn 

matter physics since their extraordinary edge and surface states pj, |2J . Following a series 
of theoretical predictions (3 6 1 , a variety of two-dimensional j3| and three-dimensional jsl- 



ri| TI materials have been realized in recent experiments. The helical spin structure of 



Dirac electrons in the gapless strong TIs acquire a spin-orbit induced nontrivial Berry phase 
of 7r after a 2tt adiabatic rotation along the Fermi surface, which results in prohibition of 



backscattering and the weak anti 
sured by transport experiments 



ocalization (WAL). The WAL effect in TIs has been mea- 



12 



la ]. Through gradually doping Cr magnetic elements 



in Bi 2 Se3 to open and increase the Dirac electrons's energy gap, most recently, the crossover 



from WAL 
prediction 



;o weak localization (WL) has been well observed [19] . consistent with theoretical 



20]. 



Meanwhile, the surface scanning tunneling microscopy (STM) measurements have been 
extensivel y c arried out to study the electronic properties and impurity scattering effects of 
TIs , |2l] . Compared to the conventional transport techniques that give the averaged 
signals, the STM has its own advantages in probing and even manipulating single-impurity 
and multiple- impurity scattering, from which the more precise quantum processes and mech- 
anisms become possible to reveal. In particular, due to the nature of the single Dirac cone, 
no complicated intervalley scattering events happen on the TI surface, therefore, in situ de- 
signing of specific impurity configurations and mapping of their scattered electronic states 
can be harnessed to image extraordinary quantum properties of TI surfaces. Recently, 
we proposed a spin-dependent Aharonov-Bohm (AB) interferometer [3], which consists of 
a spin-polarized STM tip and two identical nonmagnetic impurities, to probe the unique 
spin-texture related quantum scattering behaviors on gapless TI surfaces. Using this inter- 
ferometer, we have shown an interesting analogous WAL phenomenon reflected from $o/2 
periodic AB oscillations in the spin-resolved local density of states (LDOS), which is absent 
in the conventional metal surface systems. 

In the present paper, we extend to consider the spin-dependent AB effect on gapped TI 
surfaces. If the Dirac fermions gain mass by coupling to a magnetic exchange field, a gap will 
be induced in the system, breaking time reversal symmetry in scattering amplitudes, and as 
a consequence scattering and transport properties for gapped and gapless systems should be 
significantly different, as witnessed by the above-mentioned WAL to WL crossover. Inspired 
by this idea, thus we study the effects of finite gap on the AB oscillations in spin-resolved 
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LDOS that we previously initiated. The AB oscillations in the real-space LDOS maps are 
arisen from the magnetic flux threaded through the time-reversed self-crossing loops (see 
Fig. [T]). We find on one hand that the total LDOS exhibits the AB oscillation with a period 
of &Q=hc/e in both gapless and gapped cases. On the other hand, while the spin-resolved 
LDOS shows $o/2 periodic AB oscillations on the gapless TI surface, as we have reported 



before 



22], however, with increasing the energy gap, the $o/2 period of AB oscillations 



in spin-resolved LDOS gradually disappears and at the same time the $0 period becomes 
clear. In other words, the crossover from analogous WAL to WL on TI surface can be 
well observed using our spin-dependent AB interferometer instead of the complicated low- 
temperature transport measurement. This crossover is consistent with the evolution of the 
Berry phase with increasing the TI surface-state gap. 

The TI surface, on which two nonmagnetic impurities are adsorbed, is described by a 
low-energy effective Dirac Hamiltonian 



H = H + V(r) 



(1) 



where 

H =hv f (k x a y -k y a x ) +Aa z , (2) 

with Vf (^4~5xl0 5 m/s for Bi2Te3-family TIs as numerically used throughout this paper) 
being the Fermi velocity and 2A the energy gap of massive Dirac fermions, which is absent 
in the massless limit. 

2 

V(r)=Y,Ui<T 6(r-ri) (3) 

1=1 

denotes the potential of two impurities located at r\= (—d/2,0) and r 2 = (d/2,0) with Ui the 
potential scattering strength for scalar impurities. <j is the 2x2 unit matrix. 

The unperturbed real-space Green's function Go(r— r',u) can be obtained from the 
Fourier transformation of Gq (k,iu)=(iu-Ho)~ 1 in the fc-space, which after a straightfor- 
ward derivation is given by 

— iu 



G„(r-r» 



( 2u,r ;i+AK)fll " w 

-i(p x tr) ■ zy/S^dlH^ 



and 



\tu\ 
2^ 

3 



diag [5+fo,5-f 



(4) 



(5) 



FIG. 1: (Color online) Surface electronic interferometer, comprising a spin polarized STM tip 
at r and two impurities at r\ and r 2 separately. The interference contributions in the LDOS is 
introduced by the electrons traveling along clockwise and anticlockwise loops enclosed by the STM 
tip and two impurities. The applied magnetic field B affects this interference via the AB effect. 

for M >A, where 5± = (l±A/u), z= v /<5 +' 5 -^ anc j t r Here, 5 is the unit 

vector of p=r—r' and H^l (z) are the Hankel functions of the first kind. 

The features we discuss are expected to be seen in the change of the real-space LDOS 
owing to the influence of magnetic flux which passes through the area enclosed by the two 
scattering paths shown in Fig. [TJ This quantity can reveal the analogous WL or WAL effect 
in TI via AB oscillatory periods in LDOS. The real-space Green's function involving the 
impurities scattering is given by Dyson equation G=G +SG, with 

SG = J dr"G (r - r"; u) V (r") G (r", r'; u) . (6) 

Following the perturbation approach, Eq. (jSJ) can be expanded to any order in the impu- 
rity potential V. Our effort is concentrated on the scattering processes of surface electrons 
with the both impurities, in which the scattering paths enclose loops [23] . Therefore, taking 
all this into account, after a long algebra calculation, we have 

5G L = G (r- n) W^o (n - r 2 ) T 2 G (r 2 - r') 

+ G (r - r 2 ) W 2 G Q (r 2 - n) T X G Q (n - r') , (7) 

where the subscript L represents the loops enclosed by the scattering paths of the surface 
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electrons, and 

Ti 

Wl = cro-TiGoiri-r^TjGoiri-rj) (8) 
with T— matrices T^— ao _ v ^, , (i, j—1,2). Equation (J7J) is a general formula describing 
the two-impurity back and forth scattering of the STM-probed quasiparticles. Thus, the 
interference information during the time-reversal scattering processes are included in this 
equation. 

In the presence of a weak magnetic field, the Green's function can be semiclassically 
approximated as 

Go (r - r', u) = a «>*g ( r _ r ', u) , (9) 

where A= (—By, 0, 0) represents the vector potential. The correction of the LDOS due to 
the magnetic flux is given by 



AN L (r,u,B) = ImTr 

7T 



5G L (r,u)-5G L (r,u>) 
ANl(r,u,B) + ANi(r,co,B), (10) 



where SGl is calculated from Eq. (J7|) with Gq. It is clear that the magnetic field affects 
the LDOS via the magnetic flux, which is easy to be obtained from the integral over the 
loops (r^r 1 ^r 2 ^ ± r), §. A (I) -dl=±Bdy /2=±§. In the present setup, we focus solely on 
the (spin- resolved) LDOS at r=(x, y), which is probed by the STM tip with the same plane 
coordinates. Since the STM tip also dually participate in composing the closed trajectory 
that the Dirac electron travels, hence under a fixed magnetic field B, the AB interference 
displays the oscillations with varying the tip position along the y direction. 

The condition Z#>A/ of the semiclassical approximation in Eq. ([9]) should be satisfied, 
here B=5 T and Ef— 150 meV are chosen from which the corresponding magnetic length 
Zb~11.63 nm while the Fermi wave length A/=10.5 nm. Also, the Zeeman splitting by the 
external magnetic field is negligibly small (typically of 0.5 meV at B=5 T for Bi2Se3 film 



24j ) compared to the strong spin-orbit coupling (SOC), and thereby is neglected in this 



discussion. In a recent low temperature transport experiment [19j on magnetically doped 
Bi 2 Se3 film, an energy gap as large as 100~300 meV near the Dirac point was observed, 
which is comparative with our choise in this work. 

By choosing suitable Fermi energy Ef and energy gap parameter A, which can be con- 
trolled in experiments, the backscattering and the crossover from the WAL to WL can occur 
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on TI surface. When the Fermi energy lies in the gap, the interference signals arising from 
the contributions of ANl (r, u, B=0) are so weak that the oscillatory ellipse features as well 
as the AB effect become ambiguous, hence we only consider Ef>A. 

When the energy gap is opened, differing from the gapless case, the 

(i,j=l,2) are not proportional to a unit matrix because 



T— matrices Tj= ,/J; . 

Gq 1 (0;cu) t^Gq 2 thus the W% matrices in Eq. are no longer diagonal, resulting in 

intractable complexity in reducing Eq. (10). To get the asymptotic expression of LDOS we 
consider the lowest order in the impurity potential in Eq. (J7|), which is given by 



5Gf = U 2 G (r - ri ) G (r a - r 2 ) G (r 2 - r') 
+ U 2 G (r - r 2 ) G (r 2 - n) G (n - r') . 



(11) 



For large distances (up/hvf^l), the Hankel functions can be approximated as 
-^o/i ( z ) ~^ z \J~~^ e ^ ZT7T ^ 4: \ then the unperturbed Green's function has a simple asymptotic 
form 



G n (r - r') 



2hv 



f 



;i2) 



p 

we obtain an explicit expression of the total LDOS as follows 



where e 1 " — E^ElMl _ Substituting this equation into Eq. (TTTj) and after a tedious derivation, 



AAf> (r,w, 



sm — + sin — #2 + sin — — h 

$ / V$o / V$o 



7T$ 

sm | — - + Vi 

$ 



sm I — § 2 + sin — 20 

$ y \% 



where G=4[/ 2 (w/8vr(%) 3 ) 3/2 (l/p!p 2 d) 1/2 cos ( 



llv f 



4 / ' 



(13) 
with 



e »^i/2 = Pi/2^ x + t y) an j p l j 2 = r—ri/2- The spin-up and spin-down LDOSs are written as 



AN { I 2)Vi (r,u,B) « G5± sin 
+ 5+5 



7T$ 
7T$ 



5_i_ sin 



7T$ 



7r$ 
o 



<5 2 sin ( — =f 2(f) 



sm | ± $2 — sin =p $i 

$ / V $o 



(14) 



Equations (13) and (14) are the main analytical result of this paper. Although the total 
and spin-resolved LDOSs are now explicitly obtained, their AB oscillation periods are still 



difficult to analytically determine due to the weight coefficients 5 + and 5- that arise from 
the energy gap. However, there exist two extreme cases: (i) One is the gapless limit where 
S + =5- = l; (ii) The other one is the large-gap limit. For these two cases, one can further 
analytically simplify Eqs. (13) and (14) and get the periods of real-space AB oscillations in 
the LDOSs of TI. 

The gapless case have been discussed in our last paper j^ . Actually, when A=0 we have 
the total LDOS quantity 

Aivj 2) (r, w,B) o a [cos (2tt$/$ ) - 1] (15) 

and its two spin compents 

AiVf m (r, oj, B) ^ a [sin (27r$/$ =F 0) ± sin <j)} . (16) 

Obviously, the spatial AB oscillation period of the total LDOS is Uo=^§j for fixed B and 
fixed impurity configuration, which corresponds to a $o period in the scale of flux. 

As clearly seen from Eq. (JTBl . comparing to the total LDOS, there occurs in the spin- 
resolved LDOS additional strong SOC induced quantum interference signature. This strong 
spin interference effect deviates the real-space AB oscillations from y period when a spin- 
polarized STM tip scans along the ^-direction on the TI surface in the presence of a fixed B. 
The AB interference period along the ^-direction can be numerically determined by solving 
the zero-point equation of the factors in Eq. ( TIB]) , i.e., sin(^^=F0)± sin 0=0. We can get 



the asymptotic roots ytnj^- (nGZ) for weak B. Therefore, the AB oscillation signals for 



AN L 



(2)tA 



Bd 

occur at with a spatial period of f^=^ (i- e -; if i n the scale of flux), 



the half of ANjj. This half period could be understood as an analog of WAL effect in the 
spin-resolved LDOS. 

Turning to the other extreme case that the gap is large enough to be comparable with 
the Fermi energy Sf. In this case, Eq. (14) is simplified to be 

Aiv| 2)t (r, a, B) oc 8% [cos (27r$/$ ) - 1] 

AN { p l (r, u, B) oc <L [cos (27t$/$ ) - 1] (17) 
near the Fermi energy. Clearly, the spin- resolved LDOS AA^ 2)m turns now to display 

(2) 

a complete $o period of AB oscillations as the total LDOS AN^ does, which is totally 
different from the gapless case. In other words, the spin-resolved LDOS clearly exhibits 
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FIG. 2: (Color online) Simulations of the AB oscillations of the electronic LDOS in Bi2Te3(lll) 
surface with two nonmagnetic impurities under a magnetic field B=5 T. The left, middle, and right 
panels correspond to the total, spin-up, and spin-down LDOSs, respectively. We choose A/Ef=0 in 
(a-c), 0.3 in (d-f), and 0.6 in (g-i) with ej=150 meV. The horizonal strips in patterns are signature 
of AB effect. The green dots denote the impurity positions. The other parameters are chosen as 
Vf=4xl0 5 m/s, d=30 nm, and U=l eV. 

an analogous WL phenomenon in the large-gap limit. Except for these two extreme limits, 
to observe the AB oscillation period of spin-resolved LDOS for intermediate values of the 
energy gap, we resort to exact numerical analysis on original Eqs. (141TT01). 

From numerical calculations, we find by tuning the ratio of A/ej that the crossover 
from half-period $o/2 to period $o emerges in the AB oscillations of the spin-resolved 
LDOS. When A/ef is small, the spacial AB oscillations of the spin-resolved LDOS still 
approximately possess a half period of $o/2, i.e., the analogous WAL effect is dominant in 
spin-resolved LDOS. However, with increasing A/ef, the $o/2 period of AB effect gradually 
disappears and at the same time the $o period becomes obvious in spin-resolved LDOS. 

The typical numerical results are shown in Fig. [2] for different values of A/Ef. The 
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horizontal strips in each panel are the AB oscillation signals in the real-space LDOS. It is 
obvious that in the case of A=0, the interstrip distance is yo—^§j—56.6 nm in the total 
LDOS as shown in Fig. EJ^a), corresponding to the <3?o period of AB oscillations in total 
LDOS. Whereas, the interstrip distance becomes ?/o/2=28.3 nm in the spin-resolved LDOS 
as shown in Figs. |2{b) and Etc), corresponding to the ^ period of AB oscillations, which 
is consistent with the above analysis on Eq. ( fl6|) . From Figs. E^d-f) corresponding to 
A/£j=0.3, we can observe that the horizontal strips in spin-resolved LDOSs patterns move 
towards the horizontal strips in total LDOS, signifying a crossover from analogous WAL to 
WL in the spin-resolved LDOS. With further increasing the ratio of A/sf, the spin-resolved 
LDOS take on $ periodic AB oscillations, see the downmost panels (g-i) in Fig. El where 
A/sf is chosen to be as large as 0.6. Much smaller values of B have also been tested in 
simulations, and the calculated AB oscillation patterns of (spin-resolved) LDOS are similar 
to those shown herein. Thus that, the AB effect can be effectively studied in a semiclassical 
way in a wide range of B. 

Our findings can be understood from the view point of Berry phase. The TI sur- 
face Dirac electrons traveling along two time-reversed self-crossing loops (r^ri^r2^ ± r) 
differentiate by a Berry phase associated with spin rotation of 2tt, which is given 
by —i J Q 27r d6k (ipk\ de k \ipk} = (1—A/e)tc. Here, the eigenstates of H are expressed as 
htaH^f- (7+, Til-e i9k ) , where 7±= v /l±A/e, tan6 k =k y /k x , and e=yj (hv F kf+A 2 . 
The two clockwise and anticlockwise loops enclosed by the STM tip and two impurities in 
Fig. 1 accumulate a Berry phase of f on a gapless TI surface due to the spin-momentum 
locking, and result in the WAL effect, which are represented as AB oscillations with half- 
period of $o/2 in the spin-resolved LDOSs. While with the opening of the gap, the Berry 
phase departures from tt. The larger ratio of A/sf causes greater deviation of the Berry 
phase from n, leading to stronger WL tendency which is consistent with the observations 
from the AB interferometer proposed here. Therefore, the spin-dependent AB interferome- 
ter shown in this paper may provide a feasible approach to study the competition between 
WAL and WL by observing the spatial AB oscillation periods in the spin-resolved LDOS 
maps. 

To experimentally verify our predictions exhibited here, the spin polarized STM tech- 
nique is required, which we believe is achievable 



261 ] . Finally, we should point out that 



25 



dephasing processes have been observed in transport investigations in Bi2Se3 and Bi2Te3 



9 



films M 



19| as well as in AB-effect studies of Bi2Se3 nanowires 18). The phase coher- 



ence length of Bi2Se3 and Bi2Te3 can be as large as hundreds of nanometers, which is 
tens times of the Fermi wave length. The characteristic distance in our setup must be much 
smaller than the phase coherence length (gMcZ^), so that it is reasonable in Fig. 2 to choose 
G?=3Onm<C/0 without taking into account the dephasing processes in the above numerical 
calculations. 

In summary, we have performed a semiclassical analysis of the spin polarized STM probed 
AB oscillations in the LDOS induced by two impurities on a TI surface. With increasing 
the surface gap of TI, the crossover from analogous WAL to WL has been found in the AB 
oscillations of spin-resolved LDOS. This phenomenon may provide an important alterna- 
tive approach to testify various extraordinary quantum wavefunction properties on the TI 
surface. 
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